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ON  NON- ADIABATIC  FLAME  PROPAGATION 


By:  A.  L.  Berlad,  C.  H.  Yang,  and  R.  D.  Rowe 
CONVAIR  SCIENTIFIC  RESEARCH  LABORATORY 
San  Diego,  California,  U.  S.  A, 


I.  INTRODUCTION 

Well  defined  existence  limits  as  well  as  burning  velocities  are  general 
characteristics  of  flame  systems.  Experimentally,  it  has  been  noted  that: 

1.  Abstraction  of  heat  from  flames  reduces  the  burning  velocity 
(references  1,  2). 

2.  Sufficiently  large  heat  abstraction  from  normally  vigorous  flames 
causes  flame  extinction  or  quenching  (references  3,  4,  etc.  ). 

3.  Well  defined  fuel-oxidant  mixture  ratio  limits  exist  outside  of  which 
stable  flame  propagation  does  not  appear  to  be  possible  (reference  5). 

4.  For  systems  studied  so  far,  the  chemical  nature  of  solid,  cold 
heat  sinks  does  not  seem  to  affect  any  of  the  flame  properties  under  consid¬ 
eration  -  the  burning  velocity,  the  quenching  conditions,  or  the  concentration 
limits. 

It  thus  appears  that  flame  quenching  and  flammability  limit  phenomena, 
and  burning  velocity  variations  with  heat  abstraction  are  interrelated 
quantities  and  that  a  theoretical  treatment  of  these  phenomena  must  simul¬ 
taneously  account  for  all  these  phenomena. 
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Recent  contributions  to  the  theory  of  non-adiabatic  flame  propagation 
(references  6,  7)  have  considered  these  phenomena  for  the  case  of  laminar, 
premixed,  gaseous  flames.  This  paper  represents  another  approach  to 
a  non-adiabatic  flame  theory  which  is  consistent  with  the  previously  stated 
general  characteristics  of  flame  systems.  ^  Non-adiabatic,  one-dimensional, 
flame  propagation  through  cylinders  is  considered  for  premixed  gaseous 
fuel-oxidant  mixtures.  The  theory  is  then  extended  to  also  include  flames  in 
"a  free  radical-containing  solid  inert”  as  well  as  "a  gaseous  highly  stirred 
reactor”.  -- 

A1  though  there  are  necessary  differences  among  the  flame  models 
chosen  for  the  latter  two  cases  the  concepts  employed  for  gaseous  premixed 
flames  are  basic  to  all  three  cases. 

II.  NON-ADIABATIC  LAMINAR  FLAMES 

IN  PREMIXED  GASEOUS  MEDIA 

~  —  ~  ~  —  —  —  -  _  _  —  —  - 

IIA.  The  Model: 

In  the  case  of  adiabatic  flames  in  laminar,  premixed  gaseous  media, 
the  following  treatment  is  customary.  Final  temperature  and  the  final 
equilibrium  compositions  are  determined  by  equilibrium  adiabatic  calcu¬ 
lations.  These  hot  boundary  conditions  are  then  used,  together  with  the 
conservation  and  reaction  rate  laws  and  the  conditions  at  the  cold  boundary 
to  determine  the  eigenvalue  for  the  problem,  the  burning  velocity 
(reference  8). 
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Real  flames  are  non-adiabatic  and  the  picture  is  quite  different. 

Consider  a  non-adiabatic  reacting  medium  to  be  contained  in  an  infinitely 
long  cylinder  of  diameter,  d.  The  cylinder  axis  is  taken  to  be  coincident 
with  the  X-axis.  The  reacting  medium  is  in  a  gaseous,  homogeneously 
premixed  state.  The  following  assumptions  are  made  concerning  the 
system: 

(a)  Heat  generation  in  the  medium  occurs  as  a  result  of  a  chemical 
reaction  which  is  non-zero  wherever  reactant  concentration  and  absolute  tem¬ 
perature  are  both  non-zero. 

(b)  A  quasi-steady  flame  propagation  state  exists.  That  is,  there  exists 
a  combustion  wave  (Figure  1)  which  travels  in  the  -x  direction  with  velocity 
uq.  This  requirement  that  (duQ/d0)  =  0  implies  that  there  exists  a  finite 
distance,  xq,  between  the  maximum  flame  temperature  and  a  “cold  boundary" 

„  It 

characterizing  unburned  conditions. 

(c)  Processes  governing  the  behavior  of  the  system  include  radiative 
and  conductive  heat  transfer  to  boundaries  as  well  as  chemical  reaction 
within  the  medium.  It  is  assumed  that  there  exists  at  every  value  of  x,  an 
average  temperature  (average  over  the  tube  cross-section)  which  may  be 
used  in  characterizing  radial  and  axial  heat  transfer  and  local  reaction 
rates. 

(d)  T  =  T  >  0  for  all  x. 

w  o 

(e)  'Burned"  conditions  exist  at  x  =  oo. 


I  IB.  Mathematical  Formulation: 


The  energy  conservation  equation  may  be  written 

7 :  0) 

where  the  chemical  rate,  ($  ,  is  taken  to  be  of  the  form 

(a.) 

Then,  by  definition 

-A  f'FJ^FJt)  (?) 


Also,  for  constant  mole  flow 
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2  2 


For  the  case  Fj  (p  c)  =  p  c  ,  equation  (5)  yields 


(6) 


"to 


The  heat  loss  termlQ^  (t)  will  be  of  the  form 

of  u)  -  ^ ^  4  /^HJ 


(7) 


Examination  of  the  above  system  of  equations  shows  that  the  conservation 

equation  can  be  combined  with  one  other,  say  equation  (6),  to  yield  a  third 

order  nonlinear  differential  equation  subject  to  the  boundary  conditions 

x  =  -x  :  T  =  T  ,  c  =  c  ,  dT/dx  =  £ 
o  o  o 

(8) 

x  =  +oo:  T  =  T'  ,  c  =  0,  dT/dx  =  0 

Neither  the  conditions  at  x  =  -x  nor  those  at  x  =  +oo  by  themselves 
guarantee  the  existence  of  a  flame  in  the  region  -x-Q—  x  <  +oo.  It  is  the 
satisfaction  of  the  conservation  equation,  point  by  point,  in  conjunction  with  the 
boundary  conditions  which  allow  this  point  by  point  satisfaction  which  pre¬ 
scribes  the  characteristic  uq  value.  It  then  appears  that  the  existence  of  a 

single  value  of  u  ,  rather  than  a  band  of  u  values,  can  result  only  from  the 
o  o 

circumstance  that  all  but  one  of  these  seemingly  allowed  u^  values  are  in¬ 
consistent  with  the  conditions  of  the  problem. 


III.  FLAME  PROPERTIES  -  NON-AD  IABATIC  FLAMES 


IN  PREMIXED  GASEOUS  MEDIA 


In  order  to  examine  the  exact  behavior  of  a  given  non-adiabatic  flame,  a 
point  by  point  integration  of  equation  (1)  must  be  carried  out  for  the  case  con¬ 
sidered.  However,  it  is  instructive  to  consider  an  approximate  approach  which 
does  not  compromise  the  essential  features  of  the  problem  but  which  does 
yield  a  very  good,  general  picture  of  the  basic  structure  and  properties  of  non- 
adiabatic  flames. 


IIIA.  Approximation  Model: 


Assume  that  the  temperature  profile  on  the  interval  ~XQ—  x  ^  +co  can 
be  approximated  by  a  three  segment  continuous  temperature  profile.  The 
first  segment,  tj,  approximately  describes  the  temperature  profile  on  the 


interval  -x  ^  x  ^  0. 


(V 


Thus,  tj  is  represented  in  terms  of  two  parameters  (t  and  xq)  which  are  to 
be  so  evaluated  that  the  conservation  equation  (1)  is  satisfied.  The  second 
segment,  t^,  approximately  describes  the  temperature  on  the  interval 


t^  is  so  chosen  that  it  approximately  describes  the  temperature  in  a  region 
where  some  heat  release  occurs  even  though  the  temperature  is  declining. 
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Here  again,  equation  (1)  is  to  be  satisfied.  The  parameter,  0(  ,  is  determined 

by  3 .  Bx  is  so  chosen  that  the  third  interval,  Bx  ^  x^+oois  one  in 
o  o 

which  virtually  no  heat  release  occurs  (e.  g.  less  than  1  percent).  Thus,  the 

region  0^  x  ^  (3xq  is  a  region  of  low  heat  release  and,  in  general,  tm 

is  quite  close  to  t  .  The  linear  approximation  of  equation  (10)  is  thus 

justified.  These  conditions  then  permit  the  rewriting  of  equation  (1)  for  the 

third  interval,  where  x  3  ^  x  £  +oo, 

or 


^2-  -zpjKh)—® 

<o 


For  a  hydrocarbon-oxygen-nitrogen  flame,  the  overall  reaction  rate,  equation 
(2),  may  be  written 

to 

and,  over  a  limited  range,  approximated  by 

$  =  k2p2c2tn  (13) 

and  n  are  so  chosen  that  both  and  d$/dt  match  at  the  mid-point  of  the 

temperature  interval  of  interest.  Thus,  equation  (6)  becomes,  for 

-x  ^  x  &  0 
o 

and  for  0  —  x  —  x  3 

r\ 1 


where  for  even  n-values: 


A  similar  expression  results  for  c_, 
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where 

Thus,  equations  (16)  and  (18)  enable  us  to  determine  the  concentrations  at 

t  =  t  and  at  t  =0Ct  ,  respectively,  in  terms  of  the  parameters  x  ,  u  ,  t  . 
mm  o  o  m 

Integration  of  equation  (1),  subject  to  the  boundary  conditions  T  =  T  , 


For  the  approximation  model,  this  equation  may  be  written  as 


Now,  the  left  hand  side  of  equation  (21)  may  be  written 


|fV  mt&sm 


and  the  right  hand  side  of  equation  (21)  may  be  written 
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The  conduction  loss  term  results  from  solution  of  the  radial  heat  conduction 
equation  (reference  9).  The  radiation  loss  term  is  based  on  the  H^O  and  CO^ 
partial  pressures  and  is  evaluated  by  use  of  the  methods  and  data  of 
reference  (10).  For  the  afterburning  region,  equation  (11)  may  be  employed 
to  evaluate  the  conduction  heat  loss.  Equation  (11)  becomes 


& 

M-2 


Solution  of  equation  (25)  yields 


If  it  is  necessary  to  consider  the  radiation  loss  in  this  afterburning  region  , 
equation  (11)  becomes  nonlinear,  but  a  linear  approximation  can  be  made  by 
modification  of  G. 


The  conservation  equation,  (21),  may  now  be  evaluated. 


') 


and  for  -x  ^  x  ^  0: 
o 


It  is  now  necessary  to  determine  (3  in  such  a  way  that 

,dt2/dx||3x  =  <<Vdx>|}* 
o  o 


and 


A  value  of  (c_/c  )  0.  01  assures  a  close  approximation  to  the  true  behavior 

M  O 

of  the  flame.  Equation  (30)  then  yields 


We  now  have  a  flame  system  described  by  five  parameters 
(xq,  u  ,  tm»  0(  ,  p  )  with  four  equations,  (28),  (29),  (31),  (32).  The  fifth 
equation  for  the  determination  of  the  five  parameters  is  obtained  in  the 
following  way.  It  should  be  recalled  that  the  temperature  profile  specified 
by  equation  (1)  is  not  satisfied,  point  by  point,  by  the  approximate  temperature 
profiles  assumed  in  the  use  of  equations  (9),  (10),  (11).  Consider  the  inte¬ 
grated  square  error  of  the  conservation  equation  at  every  point.  This  error 
results  from  the  differences  between  the  actual  and  assumed  temperature 
profiles: 

where  the  t-values  to  be  employed  are,  of  course,  those  prescribed  by  the 
profiles  developed  in  this  approximate  model.  In  order  to  determine  the  u 


a 

(33. 


13 


value  which  minimizes  the  square  error  represented  by  equation  (33),  one 

may  impose  the  condition 
00 

jL  \£*<M  =  o  &+•) 

■5U0  J 

Thus,  equations  (28),  (29),  (31),  (32)  and  (34)  define  the  five  parameters 

(x  ,  u  ,  t  ,  (V  ,  0 ) . 
o  o  m  r 

niB.  Flame  Quenching  and  Flammability  Limits: 


Equations  (28)  and  (29)  represent  the  energy  balance  within  the  most 
important  portions  of  the  flame.  Energy  conservation  must  be  satisfied  if 
the  flame  is  to  exist.  Thus,  limits  of  applicability  of  either  equations  (28) 
or  (29)  define  existence  limits  for  the  flame  system  under  consideration. 

It  is  generally  considered  that  the  minimum  cq  value  which  will  allow  a 
flame  to  propagate  through  a  large  (d  =  5cm.  )  tube  at  a  pressure  of  1  atm. 
defines  a  flammability  limit;  the  minimum  tube  size  which  will  allow  the 
propagation  of  a  flame  of  fixed  cq  defines  a  quenching  distance.  Both, 

\\  •*  W  // 

flammability  limits  and  wall  quenching  distances  ,  are  predictable  from 
the  limits  of  applicability  of  the  conservation  equation.  Consider  these  limits 
of  applicability  of  equation  (29)  for  the  case  of  various  stoichiometric, 


propane -oxygen-nitrogen  flames.  Parameters  employed  in  the  calculations 

-3 

for  propane  -air  at  p  =  1  atm.  include  p  =  1.  3  x  10  ;  k-  and  n  values  were 

O  _ 

calculated  from  the  data  of  reference  (11).  and  c^  values  appropriate 
for  the  proper  t^  values  were  used.  The  c^  values  were  chosen  to  give 
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the  proper  adiabatic  flame  temperatures.  Results  of  these  calculations 
are  given  in  Table  (1). 

It  is  significant  that  extinction  conditions  are  associated  with  a  critical 
(x  /u  )  value.  One  may  consider  this  critical  value  of  (x  /u  )  to  be  a 
critical  reaction  time  for  extinction.  Further,  it  is  to  be  noted  that  the 
radiation  loss  is  relatively  insensitive  to  tube  diameter  variations,  whereas 
the  radial  conduction  loss  varies  as  d  Consequently,  the  radiation  term 
predominates  in  determining  the  flammability  limits;  the  conduction  term 
predominates  in  determining  quenching  distances.  An  important  fact  illus¬ 
trated  in  Table  (1)  is  that  critical  (x  /u  )  values  decrease  with  increasing 

o  o 

O^/N^  ratios  and  increasing  t  values.  This  explains  why  a  propane- 

oxygen  flame  quenches  at  a  much  higher  t  than  does  a  propane-air  flame. 

Further,  it  is  important  to  note  that  significant  quantities  of  reactants  are 

still  unburned  at  t  =  t  ,  particularly  under  quenching  conditions.  For  normal, 
-  m 

near -adiabatic  flame  propagation,  (C|/CQ)  i®  ®m^l. 

It  is  of  interest  to  note  the  relation  of  equation  (29)  to  previously  pro¬ 
posed  semi-empirical  quenching  equations.  Equation  (29)  reduces  to  a 
form  similar  to  those  proposed  in  references  (12)  and  (14)  if  it  is  assumed 
that,  (a)  radiation  losses  are  negligible  and  (b)  all  the  fuel  is  burned  prior 
to  t  =  t  .  For  these  simplifying  assumptions 


EXTINCTION  LIMITS  FOR  VARIOUS  STOICHIOMETRIC  FLAMES 


**>’11 1 . 


Thus,  one  should  not  expect  equation  (35)  to  present  a  completely  acceptable 
description  of  extinction  phenomena. 

JTV.  OJHER  NON-ADIABATIg .FLAM ES 

The  approach  employed  in  the  description  of  the  structure,  behavior 
and  stability  of  a  premixed  laminar  flame  may  be  extended  to  other  flame 

systems.  In  this  section,  abbreviated  discussions  of  such  approaches  are 

\\  ✓ 
presented  for  the  cases  of  Flames  in  a  Highly  Stirred  Reactor  and 

fl  // 

Flame 8  in  a  Free  Radical  Containing  Solid  Inert  . 

IVA.  Highly  Stirred  Reactor: 

VawA  — * — ‘ - 

The  combustion  behavior  of  a  highly  stirred  reactor  consisting  of  the 
volume  contained  between  concentric  spherical  shells  has  been  investigated 
experimentally  and  theoretically  (references  11,  15,  16).  The  assumption 
of  temperature  and  concentration  homogeneity  as  well  as  negligible  heat 
loss  is  generally  made  in  describing  the  reactor  behavior.  In  this  work, 
a  one  dimensional  reactor  is  considered  in  which  the  restriction  of  homo¬ 
geneity  of  temperature  and  concentration  is  relaxed.  In  addition,  the  effects 
of  heat  losses  on  the  stability  limits  are  examined. 

Consider  the  stirred  reactor  to  be  a  section  of  a  long  cylinder,  defined 
by  the  cylinder  of  temperature  Tq  and  diameter,  d,  and  two  porous  surfaces 
of  separation  x^  and  temperature  Tj.  In  the  region  0  6  x  £  x^  non-laminar 
mixing  is  taken  to  be  very  vigorous  and  of  such  a  nature  as  to  assure: 
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(a)  Volumetric  scales  of  temperature  inhomogeneities  which  are 
very  small  compared  tdjf  x^d^. 

(b)  A  time  average  temperature,  Tj  which  is  uniform  throughout 
Ot-xt  Xj. 

It  then  follows  that  there  exists  a  range  of  instantaneous  temperature  and 
concentration  distributions,  on  the  interval  0  £  x  :£  Xj,  which  will  satisfy 
the  requirements  listed  above  as  well  as  the  requirement  that  t  =  tj  and 
c  =  c  at  x  =  x..  Consider  the  average  rate  in  the  region  0  ^  x^ix., 

"  I  X  ■ 

(fl  (d&)  ~  A  A  ~  (/£q 


As  before, 

Yields 

A  - 


thus 


c  = 


Vv 


Now,  proceeding  as  in  the  previous  section,  one  obtains 


*****  ■■  #.*,*« 


The  conservation  equation  may  then  be  written 


+•  'kT(6-i)j)i1i-26/£l  +  ^r)crT 


AA0  V-  S,  J 


N 


where  ^  1,  permitting  neglect  of  reabsorption  effects. 

8 

In  general,  a  temperature  profile  similar  to  that  shown  in  Figure  (2) 

may  be  expected.  First,  consider  the  behavior  of  the  adiabatic  case: 

Neglect  of  heat  loss  processes  for  the  case  of  stoichiometric  propane  air 

reaction  at  T  =  300°K  and  p  =  1  atm.  results  in  the  allowed  values  of  u  and 
o  r  o 

(cx  /cq)  shown  by  the  solid  curves  of  Figure  (3).  The  value  kj  =  1.  455  x  10* 

is  deduced  from  the  data  of  reference  (11)  where  a  rate  of  0.  124 
3 

gm(fuel)/(cm  )(sec)  and  tj  =  467 are  used.  This  is  the  same  value  used  as  a 

starting  point  in  the  calculations  of  section  1IIB.  Interestingly  enough,  the 

3 

maximum  blow  out  velocity  of  7.  9  x  10  cm/sec  (for  a  corresponding 

to  81^>  completion  within  the  reactor)  approximately  corresponds  to  the 

same  high  temperature  reaction  time  available  in  the  reaction  zone  of  a 

laminar  flame,  where  x.-'^O.  05  cm  and  u  40  cm/sec. 

1  o 

For  the  non  adiabatic  case,  a  knowledge  of  as  well  as  the  emissive 


properties  of  the  reacting  medium  are  needed.  To  illustrate  the  effect 
of  heat  losses,  consider  the  following  case.  A.,  must  be  associated  with 
non-laminar  transport  processes  of  very  high  order;  thus  a  value  of 
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A.  t  1  x  10  cal/(cm)(sec)(  K)  is  not  unreasonable.  It  is  supposed  that 
purely  molecular  transport  prevails  outside  of  the  regime  0£r-  x  6=.x.y  The 
data  and  methods  of  reference  (10)  are  again  employed  to  calculate  the 
emissive  properties  of  the  hot  gases.  The  calculated  uq  and  c^/c^  values 
as  a  f  mction  of  tj  for  the  non-adiabatic  cases  treated  are  given  by  the  broken 
curves  jf  Figure  (3).  Of  particular  interest  is  the  existence  of  critical 
values  of  uq  below  which  this  reactor  cannot  operate.  As  expected,  velocity 
limits  exist.  In  addition,  as  the  reactor  diameter  is  reduced  from  5cm  to 
lcm,  it  becomes  apparent  that  a  quenching  limit  exists;  that  is,  there  exists  a 
critical  diameter  below  which  this  reactor  will  not  operate.  In  theory,  the 
low  temperature  branch  of  the  uq  vs.  t^  curve  may  exist  -  however,  it’ s  exist¬ 
ence  under  the  conditions  considered  has  not  been  demonstrated  experimentally. 
If  the  low  temperature  branch  is  unstable  then,  for  the  hypothetical  case  con¬ 
sidered,  quenching  occurs  for  a  tube  diameter  just  below  d  =  lcm. 


IVB.  Flames  in  Free  Radical  Containing  Solid  Inerts 

■  ■  — 1  »■  ■  ,■  "■  ■  ■  ■  —  ■  i  .  i  ■— 

Consider  a  long  cylinder  enclosing  a  free -radical  containing  solid  inert 
(for  example,  H-atoms  in  Argon).  Radical  recombination  may  occur  either 
in  the  volume  or  at  the  walls.  As  in  the  case  of  flame  propagation  in  pre¬ 
mixed  gases,  the  boundary  conditions,  the  transport  properties,  the  heat 
loss  mechanisms  and  reaction  rates  will  determine  the  flame’ s  existence, 
structure  and  behavior.  Flame  phenomena  in  such  systems  have  been  quali¬ 
tatively  observed  (reference  17).  Of  utmost  importance  to  the  flame 
behavior  of  such  a  system  is  the  form  of  the  reaction  rate  law.  There  exists 
little  information  on  the  kinetics  of  diffusion  controlled  reactions  of  the 
sort  under  consideration.  However,  the  general  form  of  the  rate  law  may 
be  deduced.  Seitz  (reference  18)  has  proposed  that  the  jump  frequency  of 
a  diffusing  atom  in  an  inert  lattice  be  given  by 

*  £$  *■)  M 
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Then,  the  number  of  jumps  per  (cm^)(sec)  is  va'  nB  .  Assuming  that 

E.»  E  ,  and  taking  the  probability  of  finding  a  similar  free  radical  at 
8  r 

an  adjacent  site  to  be  £  *  yields 


for  the  reaction  rate. 


In  view  of  the  fact  that  we  may  write  equation  (43)  in  the  form 


it  follows  that  the  above  form  for  the  reaction  rate  law  implies  no  changes 
in  the  basic  flame  model  developed  for  the  case  of  premixed  gases.  It 
may  be  necessary,  however,  to  account  for  radical- radical  wall  recombina¬ 
tion  rates,  where  the  diffusion  coefficient  has  the  form  given  by  Seitz  and 
where  the  method  of  reference  (19)  may  be  used  to  calculate  the  influence 
of  walls  on  radical  lifetimes.  Some  data  (reference  20)  exist  on  thermal 
conductivities  of  molecular  crystals.  In  addition,  the  high  energy  release 
per  reacting  pair  suggests  a  high  photon  yield  per  reaction.  This,  coupled 
with  the  fact  that  systems  of  interest  are  generally  cryogenic  implies  that 
the  dominant  heat  loss  parameter  in  the  reaction  zone  of  the  flame  may  be 
radiative. 
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Evaluation  of  these  relations  awaits  detailed  examination  of  the  existence 
limits  as  well  as  burning  velocities  associated  with  the  systems  of  interest. 

In  addition  to  the  above  considered  flame  phenomena,  there  exists  the 
usual  phenomenon  of  a  space  distributed  slow  reaction  with  the  attendant 
possibility  of  an  explosion  (reference  19). 

Thus,  for  the  case  of  flame  propagation  in  free  radical  containing 
solid  inerts  it  is  clear  that  the  boundary  conditions,  conservation  equations 
and  rate  laws,  can  be  simply  incorporated  into  the  non-adiabatic  flame 
theory  developed  for  the  case  of  flames  in  premixed  gases.  Burning  velo¬ 
city  measurements  under  varying  conditions  of  temperature  and  concentra¬ 
tion  are  then  capable  of  yielding  activation  energies  and  diffusion  coefficients 
pertinent  to  processes  involved. 

V.  CONCLUDING  REMARKS 

1.  The  above  treatment  can  be  generalized  to  take  account  of  diffusion 
processes,  changing  molar  flow  rates,  more  realistic  chemical  rate  laws, 
etc.  Simplifications  made  in  the  above  treatment  were  for  purposes  of 
convenience  and  do  not  compromise  the  essential  features  of  the  flame 
system. 

2.  One  may  derive  an  adiabatic  flame  theory  as  a  speqial  case  of  the 
non-adiabatic  theory  treated  herein. 

3.  Extinction  limits  associated  with  a  given  flame  are  characterized  by  a 
critical  reaction  time,  x  /u  . 

n*  a 


4.  Near  extinction,  significant  quantities  of  unburned  reactants  may  exist 

at  t  =  t  . 
m 

5.  The  tedious  calculation  of  the  burning  velocity  is  indicated,  but  not 
included  in  this  paper. 

6.  The  success  of  this  model  in  predicting  the  extinction  limits  for  flames 
covering  a  wide  spectrum  of  reactivities,  tube  sizes  and  pressures  argues 
for  its  validity. 

7.  This  model  can  easily  be  extended  to  treat  the  behavior,  structure  and 
stability  limits  of  a  gas  phase  flame  supported  by  the  evaporation  of 
premixed  solid  reactants. 
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VI  NOMENCLATURE 

constant 

mass  fraction  of  fuel 

initial  mass  fraction  of  fuel 

mass  fraction  of  fuel  at  x  =  0 

mass  fraction  of  fuel  at  x  =  Ox 

^  o 

specific  heat  at  constant  pressure 
cal/(gm)(deg) 

average  concentration  in  the  highly  stirred  reactor 
free  radical  species  fraction  in  solid  inert 
cylinder  diameter,  cm. 

Activation  Energy,  kcal/mole 

Atom-Atom  recombination  rate  activation  energy 

Activation  Energy  necessary  for  a  diffusing  atom  to  jump  across 
a  solid  lattice  saddle  point 

dimensionless  function  of  fuel  concentration 

dimensionless  function  of  temperature 

Geometric  factor  involved  in  channel  quenching  (reference  9) 

VoV 

characteristic  lattice  constant,  cm. 

Boltzman  Constant 
a  constant,  gm/(cm3)(sec) 
a  constant,  (2Mjk3)/(MRh3v^) 
a  constant,  gm/(cm  )(sec) 

3 

Heat  loss  function,  cal/(cm  )(sec) 


conduction  heat  loss 


axial  conduction  heat  loss  at  x  =  -x 

o 

radiation  heat  loss  (hot  gas  to  cold  walls) 
radiation  heat  loss  (hot  walls  to  surroundings) 
atomic  weight  of  inert  species 
atomic  weight  of  radical  species 

_3 

numerical  concentration  of  inert  species,  cm 
pressure,  atmospheres 
heat  of  combustion,  cal/gm 

3 

reaction  rate,  gms(of  fuel)/(cm  )(sec) 

3 

reaction  rate  in  the  solid,  gm(of  recombining  atoms)/(cm  )(sec) 


reaction  rate  order  constant 

O/tjHcxjP  0kj)  exp  f-  -  - 

A  o 

Temperature,  °K 

T/Tq,  dimensionless  temperature 


*i  j 


time  average  value  of  t,  stirred  reactor  case 
temperature  in  the  region  -xq  <  x  £  0 
temperature  in  the  region  O^x  -  PXQ 
temperature  at  x  =  0 
flow  velocity,  cm/sec 
burning  velocity,  cm/sec 


longitudinal  coordinate,  cm 
distance  defined  in  Figure  1,  cm. 
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Subscripts 

a 

o 

r 

A 


adiabatic 

initial  conditions,  except  as  noted 

radiation 

conduction 


Greek  Symbols 


d 


8 


0(t 


m 


pxo 


V 

6 


£ 

£ 


g 

w 


© 

A 

A 


e 

p 


number  of  saddle  points  (of  molar  energy  E  )  about  a  given  free 
radical  equilibrium  position 

a  post  t  temperature  at  which  the  flame  properties  are 
insensitive  to  the  existence  of  a  small  reaction  rate 
point  in  space  at  which  0(tm  occurs 
radiant  heat  loss  equation  constant 
error  function 
radiative  emmissivity  of  gas 

radiative  emmissivity  of  stirred  reactor  permeable  membranes 
at  x  =  0,  x  =  Xj 
time,  sec 

Thermal  conductivity,  cal/(cm)(°K)(sec) 

A0  ^  (t)3/4-  dt/(t-l) 

jump  frequency  of  a  diffusing  atom  in  an  inert  lattice,  sec  * 
characteristic  vibrational  frequency  of  an  interstitial  atom,  sec 
(dT/dx)  at  x  =  -xq 
gas  density,  gm/cm 


ftm.4 


•+---.-.*trrwn*Xtf  vMdtm&'i  vre***  v  ** 


tr 

*. 


% 


“l 

w2 


crystal  density,  gm/cm^ 
Stefan-Boltzmann  constant 
defined  by  equations  17 
defined  by  equation  19 
(1/2)  (H  +  ]f H2  +  (2G/dS  ) 
(1/2) (H  -  \/h2  +  (2G/d  ;  ) 
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Figure  (1) 
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FIGURE  LEGENDS 


Non-Adiabatic  Flame  Structure  (Qualitative) 


Highly-Stirred  Reactor,  One  Dimensional  Average  Temperature 
Profile 


Highly-Stirred  Reactor,  Velocity  Limits 


Theoretical  and  Experimental  Extinction  Limits  For  Various 


Stoichiometric  Flames  of  C^Hg/02/N2 


FIG.  2  HIGHLY- STIRRED  REACTOR,  ONE  DIMENSIONAL 
AVERAGE  TEMPERATURE  PROFILE 


CONDUCTION  /  RADIATION  ABSENT 
CONDUCTION  /  RADIATION  LOSSES 


FIG.  3  HIGHLY  STIRRED  REACTOR,  VELOCITY  LIMITS 


